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I. INTRODUCTION 



The QCD-inspired spin-flavor SUq quark model for 
the baryon-baryon interaction, proposed by the Kyoto- 
Niigata group, is a unified model for the complete baryon 
octet (Bg — N, A, £ and S), which has achieved very 
accurate descriptions of the nucleon-nucleon (NN) and 
hyperon-nucleon (YN) interactions [1-3]. In particular, 
the NN interaction of the most recent model fss2 [2] 
is accurate enough to compare with the modern realis- 
tic meson-exchange models. These quark-model interac- 
tions can be used for realistic calculations of few-baryon 
and few-cluster systems, once an appropriate three-body 
equation is formulated for the pairwise interactions de- 
scribed by the resonating-group method (RGM) kernels. 
The desired three-cluster equation should be able to deal 
with the non-locality and the energy-dependence intrin- 
sically involved in the quark-exchange RGM kernel. Fur- 
thermore, the quark-model description of the YN and 
hyperon-hyperon (YY) interactions in the full coupled- 
channel formalism sometimes involves a Pauli forbidden 
state at the quark level, which excludes the most com- 
pact (0s) 6 spatial configuration, resulting in the strongly 
repulsive nature of the interactions in some particular 
channels. We have recently formulated a new three- 
cluster equation which uses two-cluster RGM kernels 
explicitly [4]. This equation exactly eliminates three- 
cluster redundant components by requiring the orthog- 
onality of the total wave function to the pairwise two- 
cluster Pauli- forbidden states. The explicit energy de- 
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pendence inherent in the exchange RGM kernel is self- 
consistcntly treated. This equation is entirely equiva- 
lent to the Faddeev equation which uses a singularity-free 
T-matrix (which we call the RGM T-matrix) generated 
from the two-cluster RGM kernel. 

We first applied this formalism to three di-neutron and 
three-alpha systems, and obtained complete equivalence 
between the Faddeev calculations and variational cal- 
culations using the translationally invariant harmonic- 
oscillator (h.o.) basis [4, 5]. This formalism was also 
applied to the Faddeev calculation of the three-nucleon 
bound state [6] , employing the off-shell T-matrices which 
are derived from the non-local and energy-dependent 
RGM kernels for our quark-model NN interactions, fss2 
and FSS. The model fss2 yields the triton binding energy 
B t = 8.519 MeV in the 50 channel calculation, when the 
np interaction is employed for all the NN pairs in the 
isospin basis [7] . The effect of the charge dependence of 
the two-body NN interaction is estimated to be —0.19 
MeV for the triton binding energy [8]. This implies that 
our result is not overbinding in comparison with the em- 
pirical value, B t c * p = 8.482 MeV. If we attribute the 
difference, 0.15 MeV, to the effect of the three-nucleon 
force, it is by far smaller than the generally accepted val- 
ues, 0.5 ~ 1 MeV [9], predicted by many Faddeev calcula- 
tions employing modern realistic meson-theoretical NN 
interactions. The charge rms radii for 3 H and 3 He are 
also correctly reproduced. The non-local description of 
the short-range repulsion in the quark model is essential 
to reproduce the large binding energy and the correct 
size of the three-nucleon bound state without reducing 
the T>-state probability of the deuteron. 

Here we apply our quark-model NN and YN interac- 
tions to the hypertriton (^H) with the small separation 
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energy of the A-particle, B A cxp = 130 ± 50 keV. Since 
the A-particle is far apart from the two-nucleon subsys- 
tem, the on-shell properties of the AA~ and T,N inter- 
actions are expected to be well reflected in this system. 
In particular, this system is very useful to determine the 
relative strength of So and 3 Si interactions in our frame- 
work. We will be able to finetune the quark model in- 
teraction to the hypertriton binding energy. This enables 
firmer quark model predictions for the 1 Sq and 3 Si phase 
shifts. In fact, Ref. [10-12] showed that it is not at all 
trivial to bind the hypertriton as most meson-theoretical 
interactions fail to bind the hypertriton, except for the 
Nijmegen soft-core potentials NSC89 [13], NSC97f and 
NSC97c [14]. It is also pointed out in Rcfs. [10, 12] that 
a small admixture of the T,NN components less than 1% 
is very important for this binding. 

In the next section, the Faddeev equation of the hy- 
pertriton system, using the quark-model AW and YN 
RGM kernels, is discussed, paying a special attention to 
the AN-HN T-matrix. The results is given in the third 
section and the summary in the last section. Appendix 
gives some essential points to derive the Faddeev equa- 
tion, whose solution exactly satisfies the orthogonality 
conditions. 



II. FORMULATION 

We start from the three-cluster equation of the ANN- 
T,NN system interacting via two-cluster RGM kernels 



P[E- 



Ho - V^ieu) V£ GM (S 31 ) - v£ GM (e 32 )] 



xP* = , 



(1) 



where E is the negative three-body energy measured 
from the ANN threshold, and the free Hamiltonian 
Hq = H' + A77J3 is composed of the kinetic-energy oper- 
ators H' = hs\ + h 2 etc. and the mass term A 7713. In the 
following, the two nucleons are numbered 1 and 2, the A 
or £ is numbered 3. The equation actually implies the 
2x2 matrix form and P^S is the two-dimensional vectors 
composed of the upper component with ANN configu- 
ration and the lower one T,NN. The mass term Am3 is 
therefore a diagonal matrix whose matrix elements are 
zero for the ANN channel and AM SA = M s c 2 - M A c 2 
for the T,NN channel. Here M\ and are the A and 
£ masses, respectively. The RGM kernel V^ QM {e lj ) = 



Vb 



Gij + EijKij consists of the direct potential Vbi, , 



the sum of the exchange kinetic-energy and interaction 
kernels, Gij — Gfj + Gjj, and the exchange normalization 
kernel Kij multiplied with the center-of-mass energy Sij 
of the ij subsystem for the relative motion. These are 
also 2x2 matrices. For example, 
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The two-cluster RGM equation is expressed as 
[e 3 i-^3i-^ GM (e3i)]x-0. 



(3) 



The necessity of the projection operator P in Eq. (1) is 
related to the existence of the eigenstate of the AN-'EN 
normalization kernel; -K31IU31) = 7IU31} with the eigen- 
value 7=1. This is the most compact (0s) 6 spin-singlet 
configuration with the flavor SU3 quantum number (ll) s ; 
|u 31 ) = |(0s) 6 ; (U) 3 1 S ). We seek for the Pauli-allowed 
state of the ANN-JINN system by diagonalization 



K)(«3i|*A> = A|tf A ) 

1=1,2 



(4) 



in the antisymmetric model space, fy\(123) = — l i>\(213), 
which we denote as ty\ £ [11] in the following. The pro- 
jection operator on the Pauli-allowed space, P, is defined 
by taking the model space spanned by the eigenvectors 
with A = 0: 



P=£|*A><*> 



(5) 



A=0 



The three-cluster Faddeev equation, which is com- 
pletely equivalent to Eq. (1), is derived by using some 
nice properties satisfied by P, which is briefly discussed 
in Appendix. The total wave function of the hypertriton 
system is expressed as a superposition of two independent 
Faddeev components ip and <p: P^f = ijj + (1 — P\ 2 )<j) with 
ip £ [11]. The coupled-channel Faddeev equation reads 



yj = G (E)T 12 (E,e 12 )(l - P 12 )<P , 

<j> = G (E)f 31 (E,S 31 )(lP-P 12 <j>) , 



(6) 



where P12 is the exchange operator of particle 1 and 2, 
Gq(E) = l/(E — H ) is the three-body free Green func- 
tion for the negative energy E < 0, Ti 2 (E, £12) is the 
AW T-matrix in the three-body space, and T 3 i(E, £31) 
is the redundancy-free AN-EN T-matrix in the coupled- 
channel formalism. These T-matrices are generated 
from the NN and YN RGM kernels V^ M {e NN ) and 
Vyn M ( £ yn), respectively, according to the prescription 
essentially given in Ref. [4]. The energy dependence in- 
volved in these kernels is treated self-consistently by cal- 
culating the matrix elements of the two-cluster Hamilto- 
nian, which will be discussed in some details. 

For the NN sector, T 12 (E,e 12 ) in the ANN and 
T,NN spaces are given by the two-body T-matrix t(u>, e) 
through a simple replacement of the starting-energy ar- 
gument: 



T\ 2 (E, £12) 
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NN u 
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u NN 
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^n n — Inn [ E — AM^a — 



M N 4Ce 
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where Mm is the nuclcon mass, C/3 = Mp/M^ is the 
mass ratio for (3 = A or S, and q is the momentum of 
the residual A or E. The NN relative energy in the 
three-body space, Enn, is determined from 



e NN = (Py\h NN + V^ M (e NN )\P^) 



(8) 



which is actually the sum of the ANN and T,NN com- 
ponents. 

For the YN sector, the situation is more complicated 
since the AN-HN coupling involves a complete Pauli for- 
bidden state |u 3 i) in the two-dimensional space and the 
difference of the momentum dependent starting energies 
in the AN and UN channels is not constant because of 
the two types of reduced masses between YN and N. In 
fact, the equation satisfied by T 3 i(E, e 3 i) can never be 
reduced to the two-cluster coupled-channel equation in 
the free space. We request 

T 31 (E,e 31 ) = V* GM (e 31 ) 

+V£ GM (e 3 i)G (E)T 31 (E,e 31 ) , 
f 3 i(E,e 31 ) = T 3l {E,e 3l ) - (E - H )\u 31 ) 
x[E- e 3 i - h- 2 - Am 3 f 1 (u 31 \(E - H ) . (9) 

However, this equation is not actually satisfied, since the 
derivation is based on the unrealistic assumptions 

1) [E — £31 — hi — A7713] is channel independent, 

2) [£31 - h 31 - V^ GM (£ 3 i)] |u 3 i) = is satisfied. 

The second condition is only approximately satisfied 
since, in the strict RGM framework, it is not possible 
to use the empirical internal energies of clusters and re- 
duced masses for the relative kinetic-energy operators. 
In fact, the correct expression for 2) is 



[e - h 



Vi GM (eo)] |«3i> = 



(10) 



where £0 = £an — E[ nt with E- mt being the calculated in- 
ternal energies, and ho uses the calculated reduced mass 
unlike h 31 . We can choose (-Eint) a = f° r tne AAW 
channel, but (i?; n t)s — AM^a is only approximately sat- 
isfied in our quark models, fss2 and FSS. This difficulty 
also takes place when we try to derive the basic orthog- 
onality condition of the T 3 i(E, e 3 i) 



(u 31 \ l + G {u)T 31 {E,e 31 ) 







(11) 



which is essential to yield the orthogonality of the to- 
tal wave function through the second equation Eq. (6); 
(u 3 i|V +(l-iW> =0 (i = l, 2). 

Fortunately, these problems are completely solved by 
simply adding a small correction term to the RGM ker- 
nel, which is a procedure developed in Ref. [15] for mak- 
ing it possible to use the empirical values of the inter- 
nal energies and reduced masses in the RGM formalism. 
In Ref. [15], we have slightly modified the original two- 
cluster RGM equation [e - h - V^ GM (£ )] x = and 



considered the following RGM equation in the OCM (or- 
thogonality condition model) form: 



A[£-^ xp -y RGM (£o)]A X = 



(12) 



From here on, we omit the subscript 31 or YN as much as 
possible, in order to simplify the notations. For example, 
VrgmM = V£ GM (e ) and A = A31 = 1 - |u)(u| with 
\ u ) = \ u 3i) in Eq- (12). Furthermore, £ = e\N — E^ in 
Eq. (12) uses the empirical internal energy, E?^ — Am 3l 
and the relative kinetic-energy operator, ft.o exp — h 3 i — 
— (ft 2 /2/u exp )V 2 , uses the empirical reduced mass, £i oxp = 
Maw or Meat- Cn the other hand, we need to use £0 with 
the calculated internal energies in the RGM kernel. It is 
shown in Ref. [15] that Eq. (12) is converted to Eq. (3) 
by simply adding AG to Vrgm^o); i-e., V^ GM (£ 3 i) = 
Vrgm(£o) + AG. The explicit expression of AG is given 
in the paper (or Ae — case in Eq. (16) below). We 
use the same idea to eliminate the channel dependence 
of [E — £31 — hi — Am 3 ] in 1). Let us start from the 
off-shell extension of the AN-'EN coupled-channel RGM 
equation in the following OCM form: 

A [w - V xp - Vr GM (£o)] A X + |u)(w - e)(u\ x ) = . 

(13) 

The solution of this equation satisfies the orthogonality, 
( u \x) — f° r ^ =/= £■ The various energies in Eq. (13) 
are usually channel dependent; namely, if we use the la- 
bel (3 = A or S to specify the AN or UN channel, the 
diagonal matrix elements of these energies are given by 



£0 



= E 



C/3 + 2 q2 



2M N C/3 + 1 

h 2 /C/3 + 2 



2Mat 

Eint i 



Ca + 2 



C/3 + 1 CA + 1 



(14) 



where E is the negative three-body energy, q this time 
is the momentum of the residual nucleon. We can prove 
that Eq. (13) is equivalent to the following Schrodingcr- 
type RGM equation: 



[h c 



V, 



mod 
RGM 



(e)]x, 



(15) 



where a newly defined RGM kernel Vrgm( £ ) ^ s gi yen by 



^RGM(e) = ^ RGM (£ ) + AG, 

AG = A (AE int + Ae + Ah ) A 
- (AE int + Ae + Ah a ) , 
AE int = £f n xp - E int , Ah = h cxp - h , 
h 2 ( 1 1 

Ae 



2Mat VC/3 + 1 Ca + 1 



(16) 



The RGM T- matrix T 31 (E,e 3 i) is therefore formulated 
for this modified RGM kernel V^?^(e). By repeating 
the same process as to derive Eq. (9) with respect to 
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v rgm( £ )> we can nnd that E q- ( 9 ) ( and also E q- ( n )) is 

just valid if we replace (see Appendix for details) 



[E - e 3 i - h,2 - Am 3 ] 



[E - e AN - 



(17) 



In order to determine san or eo, we approximate w in 
Eq. (15) as 



U! = 



^ 2M JV C A + l y 





ean 

£S7V 







(18) 



with keeping the relationship s^n = £an — AMsa in the 
free space, and first calculate a AN-Y^N averaged value 



e YN = (P*\hZ P + V$3 l (e)\P*) . 



(19) 



Then, from e Y N = £aw (*|*> A + £sw (*|*) s and 
AM S a = san - £sn, we find 

£AN = £YN + AMea P E , 

esw = £yn - AM EA (1 - P s ) , (20) 

where P s = is the probability of the ~ENN com- 

ponent admixed in the hypertriton wave function. 



TABLE I: The channel quantum numbers of the hypertri- 
ton included in 15-channel Faddeev calculation with S- and 
£>-waves. For the NN-Y channels, / and j are coupled to 
the total angular- momentum | + and the isospin coupling is 
[(i§)i;t]0. For the YN-N channels, these are [iij'2]! and 
i]0. The isospin value t specifies the hyperon species; 
i.e., Y = A for t = and Y = E for t = 1. 
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0(1) 
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3 S! 


( D 2-)§ 


0(1) 


3 


3 Ol 


(sm 


0(1) 
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3 d 


(Dm 


0(1) 


5 


^0 


(sm 


0(1) 



III. RESULT 

The angular-momentum states of the NN-Y channel 
is specified by \[(Xs)I, (t\)j]\; [(\\) t, t]0), where (As)/ 
stands for the two-nucleon state, 2s+1 Aj, and t is the 
isospin value t = for A and t = 1 for S. Due to the anti- 
symmetric property of the two nucleons, (— ) A + S +* = — 1 ; 
we find that the 3 E and x O states contribute to the ANN 
channel and the X E and 3 states contribute to the T,NN 
channel. For the YN-N channel, the angular- momentum 
states are specified by \[(£isi)li, (t-2\)ji]\\ [(t\) | 5 |]0). 
Since the isospin of the hypertriton is zero, only the 
isospin T = 1/2 sector of the AN-HN interaction con- 
tributes to the hypertriton calculation. All the partial 
waves of the orbital angular-momentum are possible for 
each of the ANN and T,NN channels, which makes the 
number of channels for a particular partial-wave trunca- 
tion just three-times larger than in the triton Faddeev 
calculation. The hyperon species of the YNN chan- 
nels are uniquely specified by the isospin value t = 
or 1. For the orbital part, the parity conservation re- 
quires (— l) A+£ = (— l) £ i+ £ 2 — 1 The channel truncation 
is specified by the maximum value of the total angular 
momenta of the pairwise baryons, I and I\ , which we de- 
note J. As an example, all the channels of the standard 
15-channel calculation with S- and D-waves are listed in 
Table I. 

For the numerical calculation, we discretize the con- 
tinuous momentum variable p (or q) for the Jacobi co- 
ordinate vectors, using the Gauss-Legendre m- (or n 2 -) 
point quadrature formula, for each of the three inter- 
vals of - 1 far 1 , 1 - 3 far 1 and 3 - 6 far 1 . The 
small contribution from the intermediate integral over 
p beyond po — 6 far 1 in the two-body T-matrix cal- 
culation is also taken into account by using the Gauss- 
Legendre 713-point quadrature formula through the map- 
ping p = po + tan{7r(l + x)/A}. [These n 3 points for p 
are not included for solving the Faddeev equation Eq. 
(6), since it causes a numerical inaccuracy for the inter- 
polation.] The momentum region q = 6 far 1 - 00 is 
also discretized by the n 3 -point formula just as in the 
p discretization case. We take ni-n2-n 3 =10-10-5 as is 
used for the triton Faddeev calculation in Refs. [6, 7], for 
which well converged results are obtained. The partial- 
wave decomposition of the two-cluster RGM kernel is car- 
ried out numerically using the Gauss-Legendre 20-point 
quadrature formula. The modified spline interpolation 
technique developed in Rcf. [16] is employed for simpli- 
fying the treatment of the rearrangement of the Jacobi 
momentum coordinates. The Faddeev formalism with 
two identical particles or clusters is discussed in Ref. [17], 
together with some formulas for calculating the matrix 
elements of the two-cluster Hamiltonian. For the diago- 
nalization of the large non-symmetric matrix appearing 
in solving Faddeev equations, the Arnoldi-Lanczos algo- 
rithm developed in the ARPACK subroutine package [18] 
is very useful. 

Table II shows the results of the Faddeev calcula- 
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TABLE II: Results of the hypertriton Faddeev calculations 
by fss2 and FSS. The momentum descritization points are 
ni-n2-ri3=10-10-5 with qi-g2-<i , 3 =1-3-6 fm -1 . The calculated 
deuteron binding energy is Ed = 2.2247 MeV for fss2 and 
2.2561 MeV for FSS (e° xp = 2.2246 MeV). The heading E 
is the hypertriton energy measured from the ANN thresh- 
old; Ba is the A separation energy; enn (ean) is the AW 
(AN) expectation value determined self-consistently; and Ps 
is the EAW probability in percent. The norm of redundant 
components, A r 2 od = (u3i|P\I/) 2 (i = 1, 2), is less than 1CP 9 . 



model 



No. of 
channels 



E 
(MeV) 



(keV) 



ENN 

(MeV) 



£AAT 

(MeV) 



(%) 





6 (S) 


-2.362 


137 


-1.815 


3.548 


0.450 




15 (SD) 


—2.423 


198 


-1.762 


5.729 


0.652 




30 (J < 1) 


-2.403 


178 


-1.786 


5.664 


0.615 


fss2 


54 (J < 2) 


-2.498 


273 


-1.673 


5.974 


0.777 




78 (J < 3) 


-2.510 


285 


-1.660 


6.014 


0.800 




102 (J < 4) 


-2.513 


288 


-1.658 


6.022 


0.804 




126 (J < 5) 


-2.514 


289 


-1.657 


6.024 


0.805 




150 (J < 6) 


-2.514 


289 


-1.657 


6.024 


0.805 




6(5) 


-2.910 


653 


-1.309 


3.984 


1.022 




15 (SD) 


-2.967 


710 


-1.433 


6.171 


1.200 




30 (J < 1) 


-2.947 


691 


-1.427 


6.143 


1.191 


FSS 


54 (J < 2) 


-3.121 


865 


-1.323 


6.467 


1.348 




78 (J < 3) 


-3.128 


872 


-1.320 


6.480 


1.357 




102 (J < 4) 


-3.134 


877 


-1.317 


6.488 


1.360 




126 (J < 5) 


-3.134 


878 


-1.317 


6.488 


1.361 




150 (J < 6) 


-3.134 


878 


-1.317 


6.489 


1.361 



tions using fss2 and our previous model FSS. In the 
15-channel calculation including the S and D waves of 
the NN and YN interactions, we have already obtained 
Ba = -e d - E( A U) ~ 200 keV for fss2. The convergence 
with the extension to the higher partial waves is very 
rapid, and the total angular-momentum of the baryon 
pairs with J < 4 is good enough for 1 keV accuracy. As 
for the converged B\ values with 150-channel ANN and 
EAW configurations, we obtain B\ = 289 keV and the 
SAW component Ps = 0.80% for the fss2 prediction, 
and B A = 878 keV and P s = 1.36% for FSS. Table III 
shows the correlation between the A separation energy 
B\ and the 1 S'o and 3 Si effective range parameters of 
FSS, fss2 and NSC89 A A" interactions. Although all of 
these AN interactions reproduce the low-energy AN total 
cross section data below pa ~ 300 MeV/c within the ex- 
perimental error bars, our quark-model interactions seem 
to be more attractive than the Nijmegen soft-core poten- 
tial NSC89 [13]. The model FSS gives a large overbind- 
ing since the 1 Sq AN interaction is strongly attractive. 
The phase-shift difference of the 1 Sq and 3 Si states at 
PA ~ 200 MeV/c is S^So) - <5( 3 Si) - 29° for FSS, while 
Si 1 S ) - <5( 3 Si) ~ 7° for fss2. Since the present fss2 re- 
sult is still slightly overbound, this difference should be 
somewhat smaller in order to reproduce the correct ex- 



TABLE III: and 3 Si effective range parameters of FSS [1], 
fss2 [2, 3], and NSC89 [13] AA interactions (Ap for NSC89) 
and the A separation energies Ba of the hypertriton. The 
values for NSC89 are taken from Ref. [12]. 



model 


a s (fm) 


r s (fm) 


a t (fm) 


r t (fm) 


Pa (keV) 


FSS 


-5.41 


2.26 


-1.03 


4.20 


878 


fss2 


-2.59 


2.83 


-1.60 


3.01 


289 


NSC89 


-2.59 


2.90 


-1.38 


3.17 


143 



perimental value Pa° xp — 130 ± 50 keV. From the two 
results given by fss2 and FSS, we can extrapolate that, 
in our quark model, the desired difference is 0° ~ 2°. 

In order to make sure that this extrapolation gives 
a good estimation, we modify the K-mcson mass of the 
model fss2 from the original value, m K = 936 MeV [2], 
to 1,000 MeV, and repeat the whole calculation. It is 
known that this modification makes the 1 So AN inter- 
action less attractive and the 3 Si more attractive. We 
obtain B A = 145 keV with P s = 0.53%. The effec- 
tive range parameters of this modified fss2 interaction 
are a s = —2.15 fm, r s — 3.05 fm, and at = —1.80 fm, 
r t = 2.87 fm. The phase-shift difference is only 1.3° and 
the total cross section of the AA^ scattering increases at 
most 10 mb at p A = 100 MeV/c from 286 mb to 296 mb, 
which is still within the experimental error bars. 

It should be kept in mind that the effective range pa- 
rameters or the .S-wave phase-shift values determined in 
this way is very much model dependent, since the Ba 
value is not solely determined by these quantities. It de- 
pends on how higher partial waves contribute and also on 
the details of the AN-EN coupling of a particular model. 
A nice extrapolation shown here is based on the similar- 
ity of the models fss2 and FSS, which have a common 
framework for the quark sector and the effective meson- 
exchange potentials. 

Table II also shows that the expectation value of the 
NN Hamiltonian, Enn, determined self-consistently is 
rather close to the deuteron energy — ed, especially in 
fss2. This feature is even marked if we decompose these 
energies to the kinetic-energy and potential-energy con- 
tributions. Table IV shows this decomposition with re- 
spect to fss2, FSS and NSC89. (For this comparison, 
we use the definition of the kinetic-energy part of the 
deuteron by h d = (Xd\h N N\Xd)/(Xd\Xd}, where Xd is the 
RGM relative wave function between the neutron and the 
proton.) In fss2, the kinetic-energy of the AW subsys- 
tem is 1.88 MeV larger than that of the deuteron, which 
implies that the AW subsystem shrinks by the effect of 
the outer A-particle, in comparison with the deuteron in 
the free space. In NSC89, this difference is even smaller; 
i.e., 1.18 MeV. These results are consistent with the fact 
that the hypertriton in NSC89 is more loosely bound 
(Pa = 143 keV [12]) than in fss2 (289 keV), and the 
A-particle is very far apart from the AW cluster. The 
SAW probability in NSC89 is P s = 0.5% [10, 12]. Ta- 
ble IV also lists the kinetic-energy and potential-energy 
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TABLE IV: Decomposition of the NN and YN expectation 
values (ejvjv and £yn), the deuteron energy (— e d ) and the 
total three-body energy E to the kinetic-energy and potential- 
energy contributions. The unit is in MeV. The results for 
NSC89 are from [10]. 



lIlULLcl 


h'NN 


+ Vnn ' 


= £nn 


FCC 
r jj 


1 Q QSfi 


- 21.303 


= -1.317 


fa«9 


ly.oi u 


- 21.032 


= -1.657 


in uvoy 




- 22.25 = 


= -1.77 


111UUL1 


ft d + v d 


= — Ed (deuteron) 


r jj 


16.982 - 


- 19.238 


= -2.256 




17.495 - 


- 19.720 


= -2.225 




19.304 - 


- 21.528 


= -2.224 


model 


flYN 


+ Vyn ■- 


= £yn 


FSS 


10.036 


- 4.602 


= 5.435 


fss2 


8.071 


- 2.671 


= 5.401 


NSC89 


7.44 


- 3.54 = 


= 3.90 


model 


(Ho) + (V) 


= E 


FSS 


27.372 - 


- 30.506 


= -3.134 


fss2 


23.860 - 


- 26.374 


= -2.514 


NSC89 


23.45 - 


- 25.79 = 


= -2.34 



decompositions for the averaged YN expectation value 
Eyn and the total energy E. The kinetic energies of Eyn 
are much smaller than those of £nn, which indicates that 
the relative wave functions between the hyperon and the 
nucleon are widely spread in the configuration space. The 
comparison of the total-energy decomposition shows that 
the wave functions of fss2 and NSC89 may be very sim- 
ilar. A clear difference between fss2 and NSC89 appears 
in the roles of higher partial waves. The energy increase 
due to the higher partial waves than the S and D waves 
is 91 keV in fss2 and 168 keV in FSS, respectively. On 
the other hand, the results in Ref. [10] imply that this is 
only 20 - 30 keV in the case of NSC89. This difference 
can originate from both of the TV TV and YN interactions. 
Since the characteristics of the meson-theoretical YN in- 
teractions in higher partial waves are a priori unknown, 
more detailed analysis of the fss2 results might shed light 
on the adequacy of the quark-model baryon-baryon inter- 
actions. 



quark level. The hypertriton is an appropriate place to 
investigate the roles of the compact SU3 (ll) s compo- 
nent, which is completely Pauli- forbidden in the 1 Sq A7V- 
5W channel coupling with the isospin T = 1/2. In order 
to deal with this off-shell effect of the quark-model in- 
teraction precisely, we have formulated a new type of 
the Faddeev equation which explicitly employs the two- 
cluster RGM kernels [4, 5]. The energy-dependence of 
the RGM kernel is self-consistently treated, by calculat- 
ing the expectation values of the two-cluster Hamilto- 
nian with respect to the obtained solutions of the Fad- 
deev equation [17]. We have especially payed attention 
to how to extend the microscopic description of the A7V- 
5W coupling in the hypertriton system without spoiling 
the effect of the Pauli-principle at the quark level. The 
present study is the second application of this formal- 
ism to the few-baryon systems interacting via the quark- 
model baryon-baryon interactions, following our previous 
one to the triton system [6, 7]. The hypertriton is well 
suited to investigate the on-shell properties of the A7V 
and EAT interactions, since the hyperon is very far apart 
from the two-nucleon cluster. 

We have found that our quark-model interaction fss2 
gives a reasonable result for the hypertriton properties, 
which is rather similar to the result of the Nijmegen soft- 
core model NSC89 [13] . The A separation energy given by 
fss2 is B\ = 289 keV, which is a little too large in compar- 
ison with the experimental value _Ba° xp = 130 ± 50 keV. 
The admixture of the SAW component is P s = 0.80 %. 
Modifying the K-meson mass of fss2 from the original 
value, m K = 936 MeV, to 1,000 MeV leads to the almost 
correct A-separation energy 145 keV with Ps = 0.53%. 
Unlike the NSC89 result, the effects of higher partial 
waves up to the G wave are rather important in the 
quark-model AW and YN interactions. If we use the 
dominant 5-wave character of the AiV interaction in the 
hypertriton system, the 1 S'o AN interaction of the model 
fss2 is slightly too attractive. It is a future problem to 
investigate whether or not a reduction in the 1 So attrac- 
tion like the modification m K c 2 = 936 MeV to 1,000 MeV 
produces a favorable feature for the level spacing of the 
+ and 1 + states of and ^He systems. The fairly 
large charge symmetry breaking in these systems is also 
an important issue to understand the hyperon-nucleon 
interactions in detail. 



IV. SUMMARY 



In this study, we have carried out the Faddeev cal- 
culations, using the recent quark-model AW and YN 
interactions, FSS [1] and fss2 [2, 3]. These are realis- 
tic interactions which describe all the available AW and 
YN data, by incorporating the effective meson-exchange 
potentials at the quark level. Since these quark-model 
baryon-baryon interactions are formulated in the RGM 
framework, they are non-local, energy-dependent, and 
sometimes involve the Pauli-forbiddcn component at the 
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APPENDIX: DERIVATION OF THE 
THREE-CLUSTER FADDEEV EQUATION FOR 
THE ANNSNN SYSTEM 

In this appendix, we discuss some essential points to 
derive the Faddeev equation Eq. (6) from the three- 
cluster equation Eq. (1). In the derivation, we extensively 
use the following properties of the projection operator P 
defined in Eq. (5): 

(i) A 3l P = PA 3l = P (i = 1,2), 
where A 3i = 1 - \u 3l )(u 3l \ , 

(ii) when * e [11], V (u 3i \$) = < — > P* = # , 

(iii) when * e [11], P* = < — ► 3 |/} , 

* = |«3l/2> - |«32/l> • (A.l) 

Using the property (i), we can replace V 3 ^ GM (e 3 i) in Eq. 
(1) by v 3l (e 3l ) = A 3l V£ GM (e 3l )A 3l or 

V 3t (E,e 3l ) = {E- H ) - A 3l (E - H a )A 3l + v 3l (e 3t ) . 

(A.2) 

We further use the property (iii) for the whole equation 
and introduce the ansatz for the Faddeev components, 
Pvfr = t\) + (1 — P\2)(j), to derive a pair of equations 



[E-H - 
[E-H 



^ GM ( £ 12)]V 



^ GM ( £ 12)(1-Pl2) 

V 3 i(P,£ 3 i)]0 = V 3 i(S,e 3 i)(^-Pi2< 



+ IW31/2) 

In the second equation, we note that 

E — H — V 3 i(E, £31) 

= A 31 [E-Ho-V* GM (s 3i ) 



A 



31 



(A.3) 



(A.4) 



and introduce the projected two-body Green function in 
the three-body space, G\ 31 (E,e 3 i), which satisfies 

G Asi (E,e 31 ) A 31 [E - H - v 31 (e 31 )]A 31 = A 31 . 

(A.5) 

This can be easily constructed through 

G\ 31 (E, e 3 i) = G V31 (E, £31) - G V31 (E, £ 3 i)|u 3 i) 
x 7 \n — n? I \\ r(«3i|G«3i(^.£3i) , (A.6) 

{U3\\G V31 {E,£ 31 )\u 3 i) 

by using the two-body Green function G V31 (E, £31) = 
[E — Ho — v 3 i(e 3 i) + iO]^ 1 in the three-body space. The 
essential equation we need for deriving the full Green 
function G 31 (E, e 31 ) = [E — H - V£ GM (£ 31 ) + iO] ^ is 
the decomposition 



(£31) = (E — e 3 i - h- 2 



Am 3 ) 



P. II,, lf ; 

-A31 (E - £31 -fi2- Am 3 ) A31 

+A 31 [P-i? -^ GM (£ 31 )] A31 

= |u 3 i) (E - e 3 i - /15 - Am 3 ) (u 3 i| 

+A 31 - ifo - ^ GM (£3i)] A31 , (A.7) 
but the last equality is actually not satisfied since 
(E — £31 — /13 — Am 3 ) is channel dependent. This dif- 
ficulty is avoided by using V^^e) in Eq. (16), in place 
of y 3 ^ GM (£ 3 i). In fact, we find that 

(E-eo-h- 2 - Am 3 ) - {h 31 - h ) + (AE int + Ae 
+Ah ) = E - e AN - (h- 2 ) A (A.8) 

is channel independent. Here £0 and (h 3 \ — ho) term ap- 
pear since |u 3 i) actually satisfies Eq. (10) and not Eq. (3). 
This makes it possible to derive our basic relationship 



G (E)T 31 (E,e) = G 31 (E,e)V^(e) = G A31 (E,e)V 31 (E,e) - |« 3 i)<«3i| + l«3i> 



1 



E - e AN - (h 2 ) j 



-(u 31 \(E-H ) 



G (E)T 31 (E, e) + \u 31 ) 



E - EAN - (h 2 ) A 

r 



(u 31 \(E-H ) 



(A.9) 



where all the kernels are defined by using V^?^(e). From tion Eq. (6) and the orthogonality condition Eq. (11). 
Eq. (A.9) we can easily prove the second Faddeev equa- 



[1] Y. Fujiwara, C. Nakamoto, and Y. Suzuki, Phys. Rev. 

Lett. 76, 2242 (1996); Phys. Rev. C 54, 2180 (1996). 
[2] Y. Fujiwara, T. Fujita, M. Kohno, C. Nakamoto, and Y. 

Suzuki, Phys. Rev. C 65, 014002 (2002). 
[3] Y. Fujiwara, M. Kohno, C. Nakamoto, and Y. Suzuki, 

Phys. Rev. C 64, 054001 (2001). 
[4] Y. Fujiwara, H. Nemura, Y. Suzuki, K. Miyagawa, and 



M. Kohno, Prog. Theor. Phys. 107, 745 (2002). 
[5] Y. Fujiwara, Y. Suzuki, K. Miyagawa, M. Kohno, and H. 

Nemura, Prog. Theor. Phys. 107, 993 (2002). 
[6] Y. Fujiwara, K. Miyagawa, M. Kohno, Y. Suzuki, and H. 

Nemura, Phys. Rev. C 66, 021001(R) (2002). 
[7] Y. Fujiwara, K. Miyagawa, Y. Suzuki, M. Kohno, and H. 

Nemura, Nucl. Phys. A721, 983c (2003). 



8 



[8] R. Machleidt, Adv. Nucl. Part. Phys. 19, 189 (1989). 
[9] A. Nogga, H. Kamada, and W. Glockle, Phys. Rev. Lett. 
85, 944 (2000). 

[10] K. Miyagawa, H. Kamada, W. Glockle, and V. Stoks, 

Phys. Rev. C 51, 2905 (1995). 
[11] K. Miyagawa, H. Kamada, W. Glockle, H. Yamamura, T. 

Mart, and C. Bennhold, Few-Body Systems Suppl. 12, 

324 (2000). 

[12] A. Nogga, H. Kamada, and W. Glockle, Phys. Rev. Lett. 

88, 172501 (2002); A. Nogga, Ph. D. thesis, University 

of Bochum (2001). 
[13] P. M. M. Maessen, Th. A. Rijken, and J. J. de Swart, 

Phys. Rev. C 40, 2226 (1989). 



[14] Th. A. Rijken, V. G. J. Stoks, and Y. Yamamoto, Phys. 

Rev. C 59, 21 (1999). 
[15] Y. Fujiwara, M. Kohno, C. Nakamoto, and Y. Suzuki, 

Prog. Theor. Phys. 104, 1025 (2000). 
[16] W. Glockle, G. Hasberg, and A.R. Neghabian, Z. Phys. 

A 305, 217 (1982). 
[17] Y. Fujiwara, K. Miyagawa, M. Kohno, Y. Suzuki, 

D. Baye, and J.-M. Sparenberg, KUNS-1910, nucl- 

th/0404071, submitted to Phys. Rev. C. 
[18] See ARPACK homepage, http://www.caam.rice.edu/ 

software/ARPACK/ 



